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1, INTRODUCTION 

Time synchronization is a crucial issue to be addressed in the implementation of WSANSs. 
Through time synchronization, all the sensors and actuators in the network will follow the same time 
reference. With the help of same time reference, WSANSs can sequence the events and model the physical 
environment sensed by the sensor nodes accurately [1]. Besides, a same time reference among the sensors is 
also influential in identify the casual relationships between events in the physical world, to eliminate 
redundant sensor data, and to generally facilitate sensor and actuator network operation [2]. 

There are many synchronization protocols available in the literature, for example, Reference 
Broadcast Synchronization (RBS) [3], Time-sync Protocol for Sensor Network (TPSN) [4], Reachback 
Firefly Algorithm (RFA) [5], Average TimeSync (ATS) [6], Time Synchronization using Maximum and 
Average consensus protocol (TSMA) [7] and Proportional-Integral clock Synchronization (PISync) [8]. 

As pointed out in [8], most of the protocols in the literature, especially those based on least-squares 
regression, have a drawback where the effect of various error sources will appear as multiplicative noise in 
the time synchronization error dynamic, which in turn cause the global synchronization error to increase 
exponentially with the diameter of the network, approximately. This drawback will then cause the poor 
performance in term of scaling properties. This flaw can be solved by using control theory based approach, 
where the synchronization is accomplished using linear feedback on the measured local synchronization 
error. Through this control-theoretic approach, the error sources will appear as additive noise in the error 
dynamic and the global synchronization error will proximately increase as the square root of the diameter of 
the network [8]. 
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The time synchronization protocol presented in this paper is a control-theoretic approach. It is called 
Time Synchronization using Distributed Observer algorithm with Sliding mode control element (TSDOS), 
which applies a distributed observer/estimator algorithm and sliding mode control element on relative 
comparative error to estimate the skew rate and relative skew rate to accomplish the purpose of time 
synchronization. 

TSDOS has the following characteristics: distributed, asynchronous, robust and light in term of 
computation effort. The rest of this paper is organized as follows. Section 2 describes the development of 
TSDOS algorithm and the convergence analysis of the TSDOS algorithm, followed by the brief details of the 
simulation setup in the last sub-section of Section 2 while the simulation results are displayed and discussed 
in Section 3. The simulation results of TSDOS are compared with another two time synchronization 
protocols from literature: ATS [6], and PISync [8] Finally, Section 4 concludes the findings. 


2. TSDOS METHODOLOGY 
2.1. Clock Modelling 

Every node t in a WSAN has its own local clock which is a monotonically non-decreasing function 
of t. Each local clock dynamic can be modeled as first order dynamics, as shown in (1) 


T,(t) = a;t + B; (1) 


where T; 1s the local clock reading, @; is the local clock skew which determines the clock speed, and f; is the 
local clock offset which is the difference between time reported by clock T; and the real time at the initial 
instant Co. 

Figure 1 displays relative clock dynamics between a pair of nodes with different skew rate and 
offset [6]. Different nodes have different skew rate and offset due to the low-cost oscillators in the nodes. 





Figure 1. The clock dynamics of two nodes, t and j, with their respective local clock skew rate, a, and clock 
offset, B. 


To achieve absolute synchronization in a sensor and actuator network, all nodes 1 = 1,2,...V must 
precisely compensate for their clock parameters a; and f; so that all clocks have same clock skew and zero 
offset error [9]. Without information on the real world absolute time t available to the nodes in WSAN, 
computing for the clock parameter a; and fh; is impossible. As reported in [10], to avoid this issue, firstly, 
we can rearrange (1) so that t is expressed as a shared entity, 1.e. 


t= Ti-Bi (2) 


Qj 


and then substitute (2) into (1) for node j so that we will reach 
Qa; Qa; 


As can be observed, as shows in (3), local clock dynamic for node j does not contain the absolute 


reference time t. The paradigm adopted in [6] is that, we may eventually obtain the clock skew rate and 
offset of all the nodes in a network by estimating the relative magnitude = and the offset imposed 


L 
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a; Fs Fs . . A 
(B; = a B;). If these magnitudes can be estimated eventually, the clock dynamics estimates T; can be 
U 


furnished indirectly. Therefore, we may express the clock dynamics estimates of each node in terms of these 
estimated magnitudes, 


T; = Q;T; + 0; (4) 


where @; is the relative clock skew rate estimate and 6; is the clock offset estimate. As shows in (4) that each 
node has its own self-time keeping mechanism by estimating its own virtual time, T; before the estimates is 
passed to its connected neighbor node j. The final objective is to have a virtual reference 
clock, T,., 


T,(t) = a,t + 0, (5) 
to which all the nodes will refer to and finally be synchronized, 1.e. 
litizact, = t.5 Vi (6) 


The end objective and motivation is like the TSMA algorithm [7] whereby, as shows in (6) describes 
the consensus in achieving the common reference time except that TSMA adopts the maximum averaging of 
all the virtual time estimates. On the other hand, ATS algorithm in [6] adopts a kind of low pass filter 
structure from which, in this paper, we have enhanced the algorithm by imposing a sliding-mode like term 
for fast convergence in relative skew estimates and skew rate estimates. 


2.2. TSDOS Algorithm 

TSDOS algorithm is an extended algorithm from the work of ATS [6] through an enhancement by 
the sliding-mode control term. Two assumptions are deduced for the algorithm to work. First, the 
communication created between neighboring nodes is a single-hop communication which occurs at a very 
short finite time interval. Second, the communication occurs between a node pair is instantaneous. 

TSDOS consists of three main parts: the relative skew rate estimation, the skew rate estimation, and 
the offset estimation. In this work, each node I is assumed to periodically transmit a packet to all its 
neighbors with a synchronization period equal to T, i.e., the transmission instant t} 


i €T—B; 
2 a (7) 
2.2.1 Relative Skew Estimation 
Every node t tries to estimate the relative skew rate with respect to their neighbor j. The relative 
skew rate estimates in [6] is enhanced by incorporate a sliding-mode like term at the relative skew rate. 
We first define a sliding plane e, 


aa Tj (t2)—T;(€1) 
a Tj (t2)-Tj (C1) 


(8) 


from which i.e. 7; ;(t,) — 7;,;(t,) taken at time t, and time t, from the data packet for node j is compared 
with node t. Then the relative comparative error is fed into a low-pass filter structure to estimate the relative 
skew rate, 


S(t + 1) = peéi(t) + (1 — pg)sign(e,) (9) 
where p; € {0,1} is the tuning parameter similar to the effects of the filter poles coefficient, i.e. increasing 


the coefficient may increase the convergence of the filtered signal but may reduce the effects of filtering. 
A signum function forces the comparative clock pairs to the plane defined in (8). 


Theorem 1. Consider the relative skew update (9) where 0 < p; <1, the transmission events th are 
generated according to the (7), and each T; evolves according to (1). Then 


lim;-,c0+ E, (t) = 1 
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in finite-time for any initial condition €;(0). 
Proof. See Appendix A. 


2.2.2 Skew Estimation 
As stated previously in (4), a node updates its own perceived virtual time. The skew rate estimate is 
also enhanced by the sliding mode term in the filter, 


a(t + 1) = p,@,(t) + 1 — p,)éisign(4;) (10) 


where p, € {0,1} is the tuning parameter and @; is the virtual clock skew estimate of the connected neighbor 
node J. 


Theorem 2. Consider the skew update equation given by (10) with initial condition @;(0) > 0 and 0 < p, < 
1, where €;(t) is updated according to (9) and t} are defined in (7). Then 


a, = lim;_,.+ @ (t) = 1, Vi,a;(0) > 0 
in finite-time. 
Proof. See Appendix B. 


2.2.3 Offset Estimation 
Last but not least, after all the nodes are running at the same speed, the offset 6, of the virtual clock 
can be estimated by, 


6;(t + 1) = 6,(t) + 1 — p,)(4; — fi) (11) 


where p, € {0,1} is the tuning parameter and 7;,t; are the estimated virtual clock for node i and their 
neighbor node j, respectively. 


Preposition 1 (Theorem 6, [6]). Consider the offset update equation given by (9) with initial condition 
6,(0) = 0 and 0 < p, < 1,where 7;, and t5, are defined in (4) and (7), respectively. Then 


lim, cot T; (t) = 7)(0), VL E Nj 
exponentially fast. 


Finally, the virtual clock estimator, T, for each node i is then obtained by combining the skew 
estimator and offset estimator value from (10) and (11). 


2.3. Simulation Setup 

An experimental simulation is designed in this paper to verify the proposed time synchronization 
protocol, besides comparing the performance with another two control-theoretic protocols from literature, 
ATS [6] and PlISync [8]. The experimental simulations are simulated using MATLABTM SimulinkTM 
where ad-hoc deployment of 5 nodes in a WSAN is simulated. Figure 2 shows the connectivity among the 5 
nodes in the simulation. The final objective of the simulation is all the 5 nodes will eventually be 
synchronized to a same time reference. 


Figure 2. Connectivity of the 5 nodes 
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The synchronization period in this simulation is set to be | milliseconds and is running for a total 
of 100 seconds. In the simulation, cumulative Integral Absolute Error or cLAE which can be expressed as 


cIAE = YN, lejldt (12) 


Was computed for the whole nodes in the network as a measure of the computational effort (control 
effort) of the protocol. N’ denotes the total number of nodes in a network whose perceived virtual time are to 
be synchronized. Such computed index also gives an indication of the energy consumption of each node to 
achieve time synchronization, thus the lesser value is desirable. This is one of the important traits for a 
wireless sensor node as they need to conserve energy. 


3. RESULTS AND DISCUSSION 

Figure 3 illustrates the plots of virtual time estimator, T, while Figure 4 displays the plots of virtual 
time estimation error, T, for the three algorithms, TSDOS, ATS, and PISync. In general, it can be observed 
that all the three algorithms are able to synchronize the nodes. As observed in Figure 3, the virtual time 
estimate of all three algorithms are increasing linearly, which fulfill the first order dynamic stated in (5). 
In addition, from the plots of virtual time estimation error in Figure 4, it can be deduced that the error of all 
the 3 algorithms are remain constant and relatively small throughout the synchronization period. 

However, more precisely, from Figure 3 (c), it is noticeable that the line width of the virtual time 
estimate plot of PISync is quite thick compared to another two figures, Figure 3 (a) and Figure 3 (b). 
Next, from the error plot in Figure 4, 1t explained the reason of the thickness, where the error of TSDOS and 
ATS 1s converged to O (Figure 4 (a), 4 (b)), while the error of PISync is converged to a bigger range of value, 
between -6 to 12 (Figure 4 (c)), which shows the performance of PISync is unfavorable compared to TSDOS 
and ATS. 
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Figure 3. Plots of virtual time estimate of (a) TSDOS, (b) ATS, and (c) PISync 
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Figure 4. Plots of virtual time estimation error of (a) TSDOS, (b) ATS, and (c) PISync 


The cIAE of the three algorithms are summarized in Figure 5 to have a clearer comparison on the 
three algorithms. From Figure 5, it can be clearly seen that TSDOS and ATS once again outperformed 
PISync, where TSDOS and ATS have a much lower cIAE, 1.1089 x 10% and 1.6653 x 10%, respectively, 
compared to the cIAE of PISync, 6.1558 x 107. In addition, the bar charts in Figure 5 also clearly shows that 
TSDOS has the best performance which can be easily observed from Figure 5 that the bar chart for TSDOS is 
the lowest among the 3 bar charts. This is because TSDOS has the lowest value in term of cIAE, compared to 
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the other 2 protocols. As mentioned previously, cIAE also give an indication on the energy consumption of 
the protocol. Thus, we can conclude that TSDOS has reduced 33.43% and 99.99% of the energy consumption 
if compared to ATS and PISync, respectively. 
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Figure 5. cCIAE of TSDOS, ATS, and PISync 


4. CONCLUSION 

In conclusion, this paper has successfully proposed an enhanced version of distributed time 
synchronization protocol. Firstly, this protocol is claimed to be completely asynchronous because it does not 
require synchronous packets transmission between nodes. Secondly, it is a distributed protocol and has low 
computational energy requirement because each node is run the same light observer/estimator algorithm to 
estimate the virtual reference time. The proposed TSDOS algorithm has the best performance among the 
three compared protocols, where TSDOS has the lowest cIAE. Future work will involve the hardware 
implementation of the proposed time synchronization algorithm on a group of interconnected sensor 
node network. 
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APPENDIX A 
Proof of Theorem 1 
By writing the (9) in recursive formula, it follows that 
£-1 


E(t) = pzé(0) + » p? (1 — pz) sign(e;) 
h=0 
= pzé(0) + (1 — pz)sign(e,) 


(t-B;) 


Tj 


1-of 
where f = | |;and > pz = — (from finite geometric series). Since 0 < pz; < 1, then 
“Ee 


Jim, ©) = Jim, p§§ (0) + sign(e,)(1 — p#) 


= sign(é;) 
Given the properties of signum function are 


1, A>O 
sign(A) =40, A=0 
—1, A<0 


w(t) —4)(t) 


ti(t2)—t(41) 


Since €; = and T is a monotone increasing function, so €; > 0. Therefore, 


jim, &(t)=1 


and the convergence is in finite-time. 


APPENDIX B 
Proof of Theorem 2 
Following the similar notion as in Theorem 2 and Appendix A. 
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f-1 
@(t) = pi@,(0) + >» pe (1 — p,)é(t)sign(@;) 
h=0 


= p£a@,(0) + sign(@,)é,(t)(1 — p) 


From Theorem 2, lim, ,.o+ €;(t) = 1. 


a, = lim, @;(€) 


jim, pr@;(0) + sign(a)(1 — ps) 


sign(a,) 
=1 


and the convergence is in finite-time. 
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